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Abstract 

We study the time behavior of wave functions involved in tunneUng through a 
smooth potential barrier in one dimension in the semiclassical limit. We determine 
the leading order component of the wave function that tunnels. It is exponentially 
small in 1/h. For a wide variety of incoming wave packets, the leading order tunneling 
component is Gaussian for sufficiently small h. We prove this for both the large time 
asymptotics and for moderately large values of the time variable. 
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1 Introduction 



The goal of this paper is to study the semiclassical hmit of solutions to the one-dimensional 
time-dependent Schrodinger equation that involve tunneling through a simple potential bar- 
rier. Numerical simulations that illustrate our results are presented in [5]. A related wave 
packet "spawning" algorithm is also presented there. 
Specifically, we consider solutions to 

^h— = (-y^ + ^(^)j^' n;t,h)eL\R), (l.l) 

for small values of h, where the potential V is an analytic function that represents a barrier. 
Our goal is to present formulas for the part of the wave function that has tunneled through 
the barrier. 

1.1 A Qualitative Synopsis of our Results 

Our results, stated in Theorems 12.41 and 12.51 and Corollary 12.61 are quite technical, so we 
begin with an informal, qualitative discussion of a special case. 

Suppose V{x) is a very simple, analytic bump function that tends sufficiently rapidly to 
zero as x tends to +oo and — oo. For a wave coming in from the left, we choose generalized 
eigenf unctions that satisfy 

{^tkx/h _^ j^^ ^-ikx/h ^ 

T{k, h) e^'^^/^ as X — )■ +00, 

with > 0. 

We take superpositions of these, with energies below the top of the barrier V , to form 
wave packets and let them evolve. For a wide class of such superpositions, we have the 
following: 

1. If the average momentum of the incoming wave packet is ?7, then the transition 
probability for tunneling is strictly greater than |T(?7, h)^^. 

2. The average momentum of the tunneling wave packet is strictly greater than rj. 

3. The leading term for the tunneling wave packet for small h is a. complex Gaussian. 
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Our main results underpin the numerical spawning algorithm of [5] that describes semi- 
classical tunneling. The qualitative item 3 is crucial to the algorithm, whose quantitative 
information is determined numerically. 

Intuitively, items 1, 2, and 3 are easy to understand. When h is small, the function 
\T{k, h)\ increases extremely rapidly as k increases. For k near t], it typically grows like 
C exp(Q;(A; — 'r])/h), with a > 0. Thus, higher momentum components of the wave packet 
are much more likely to tunnel than the average momentum components. Items 1 and 2 are 
consequences of this observation. 

One can understand item 3 and learn more about item 2 by examining the transmitted 
wave packet in momentum space after tunneling has occurred. For example, if the incoming 
wave packet in momentum space is chosen to be asymptotic to one of the semiclassical wave 
packets (pj of 

^-uky(2n) 2-.72 (^.,)-i/2 ^-1/4 ^-1/4 H^^{k-7])/h'/^) exp{-{k-r]f/{2h)), 
then the transmitted wave packet behaves like 

C exp{a{k-r])/h) e-''^'/'-^''^ 2-^/^ (jO"'^' tt-^/^ h,-^/^ H,{{k-r,) /h}'^) ex^{-{k-r,f / {2h)). 
This equals 

^-itk^/(2n) ^ gW/(2ro 2-^72 (^■,)-i/2 ^-i/4 ^-i/4 ^^.(^^ _ r])/h^''') exp(-(A; - - af/{2 h)), 

The Gaussian factor is large only near k = rj + a, and near there, the Hermite polynomial 
behaves like its leading term, 2^ /fV/"^. This product is asymptotically another Gaussian 
with momentum near 7] + a for sufficiently small h. 

We note that the C in these expressions behaves like e~^^'^', so none of this can be 
determined by a perturbation expansion in powers of h. 

To obtain quantitative results, we must insert energy cut offs and deal with many other 
technicalities, but the description above gives an intuitive summary of our results. 

The precise statements of our results are presented in Section 2. Theorem 12.41 presents 
the very large time behavior of the tunneling wave function. Theorem 12.51 and Corollary 12.61 
describe the behavior of the tunneling wave function for all times shortly after the tunneling 
has occurred. 
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2 A More Precise Description of the Problem 



We consider wave packets that have their energy locahzed in an interval A = [Ei, E2], and 
we assume the potential V satisfies the following hypotheses: 

i) XI—)- V{x) is real analytic in the strip Sa = {z : \lmz\ < a}, for some a > 0. 

ii) There exist V{±oo) E (—00, £'1), u > 1/2, and c < 00 such that 

limsup \\Rez\'^+'' sup \V{z) - V{±oo)\ ] < c. (2.1) 

iii) For any E E A, the function V{x)—E has exactly two simple zeros, xo{E) < xi{E) 
with Xi{E) < 1. 

Under these hypotheses, we can decompose our solutions as superpositions of generalized 
eigenvectors of 

whose energy lies in A. 

Since we are interested in the tunneling process, we assume that in the distant past, the 
wave packet was a scattering state, coming in from the left of the potential barrier. Our goal 
is to describe the leading order behavior of the wave packet on the right of the potential 
barrier for small h and sufficiently large positive times. This tunneling wave is well-known 
to have exponentially small norm in 1/h. We determine its leading order component and 
show that for a wide variety of incoming states, it is a Gaussian. 

Analogous results for exponentially small reflected waves when the energy is strictly above 
a potential bump are presented in [2]. Similar results for non-adiabatic transitions in the 
Born-Oppenheimer approximation are presented in |2] and [TT] . 

2.1 Generalized Eigenfunction Expansions 

For any flxed energy E < max V{x), we let ip{x, E, K) be the solution to the stationary 
Schrodinger equation 

- h'-^ = 2iE-Vix))^ ^ p\x,E)^, (2.2) 
that we construct below. Here, 

p{x, E) = ^2{E - V{x)) > 0, for |x| > 1, 



is the classical momentum at energy E. The turning points xo{E) < xi{E), given by 
the two solutions of p{x, E) = are branch points of p{x, E), viewed as a multi- valued 
analytic function of x. We use the multivaluedness of this function in the analysis below. 
We consider wave functions 

^(x, t,h)= [ Q{E, h) e"'*^/'^ ^(x, E, h) dE, (2.3) 

J A 

for some sufficiently regular energy density Q{E, K) defined on A. Such a function is a 
solution of (II. ip under the hypotheses we impose below on the energy density Q. 

We shall derive a space-time description of the exponentially tunneling wave to leading 
order as /i — ?► for large positive times, when this wave is far enough from the potential 
bump. 

The energy densities Q{E^ K) we choose are sharply peaked at a specific value 
Eq G (-El, E2). Specifically, we consider 

Q{E, h) = e-^(^)/'"' e-'^(^)/'"' P{E, h), (2.4) 

where: 

(CI) The real-valued function G > is in C^(A), is independent of and has a unique 
non-degenerate minimum value of at Eq E {Ei, E2). This implies that 

G{E) = g{E-E^f/2 + O {{E - E^f) , where (? > 0. 
(C2) The real-valued function J is in C^(A). 

(C3) The complex- valued function P{E, h) is in C^(A) and satisfies 



sup 

e>0 



Qn 

P{E, h) 



dE' 



< Cni for n = 0, 1. 



2.2 The Specific Generalized Eigenfunction E", h) 

We first recast the eigenvalue problem for H{h) as a first order system of linear equations. 
We then use the analyticity of the potential and consider the extensions of these equations 
to the complex x-plane to perform our asymptotic analysis. The function p(-, E) has branch 
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points at xo{E) and xi{E) which we have to select in a consistent way. Note the turning 
points Xq{E) < Xi{E) are separated for all E E A. 
We initially confine x to the real axis and define 



Pr{x,E) = y/2\E-V{x)\ > 0. (2.5) 

We begin our analysis for x > xi{E), where pr{x, E) is the classical momentum. 

Suppose C satisfies the ordinary differential equation (12. 2p . For x > Xi{E), we define 

where ' denotes the derivative with respect to x. We note that $ satisfies 

th^'ix) = (p2(^%) J) H^) = Aix,E)^x). 

In order to apply the WKB method, we expand $(x) in terms of the instantaneous eigen- 
vectors of the generator of A{x, E). Specifically, we choose 

^i(x) = I I and ^,{x) = ( v ^^"'^'^^ I (2.7) 

\ ^pn{x, E) j \ - ^pr{x, E) I 

to be the eigenvectors associated with the eigenvalues p_r(x, E) and —Pr{x, E), respectively. 
We decompose $( clS 

$(x) = ci(x) e-^^-'i(^)^'«(^'^)^^/'^ y,i(x) + C2(x) e'/:i(i=)P«fe.^)'^?^A ^^^^^^ (2.8) 
where ci(x) and C2(x) are complex- valued coefficients that satisfy 

cUx) \ _ (., ^) / ^i:..^^iy^^^^yf^ \ / ,^(,) X 

4(X) y 2p^(a;, E) 1 ^-2if:^^^^pa{y,E)dy/h Q M C2(X) y 

These coefficients depend on E, and on sets where p{x, E) is real-valued, they satisfy [12] 
• |ci(a;)p — |c2(x)p is independent of x, and 



Ci(x) \ / ci(x) , 

is a solution to (12.90 if and only if is a solution to (12. 9p . 

C2(a;) / \ C2{x) 
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This decomposition allows us to write the solution to (12 ■2p for large positive x as 

C(x, E, h) = ci{x, E, h) e-^/^'i(^^)P«(2''^)'='^/'^ <^W(a;, E, h) 

+ C2(x, ;i) e'^^i(^)^«^^''')''^/'>^'^(x, /i), (2.10) 

where the E and h dependence is explicit, and (p'^p denotes the first component of (pj. 

When X is smaller than xo(-E'), so that the classical momentum again equals pji{x,E), a 
similar expansion is valid. However, we take a different phase convention when x < xo{E): 

C(x, E, h) = di(x, E, h) e-'^^oiE)Pi^^y^^^'^y/^ 

+ d2{x, E, h) e'/4(s)P«(^'^)'^^/'^ <^W(a;, E, K). (2.11) 

where the coefficients di and 62 satisfy a differential equation similar to (12. 9p . 

We need to connect the coefficients ci and C2 to the coefficients di and ^2- One commonly 
used technique consists of solving a similar equation for xo{E) < x < xi{E) and matching 
the solutions to those we just described. Instead we use the complex WKB method which 
allows us to work with just one equation, but requires an analytic framework. See [1] for 
various possible approaches. 

We only want to have an outgoing wave on the right and are not currently worrying 
about normalization, so we consider the asymptotic conditions 

ci(+oo, E, h) = 0, and C2(+oo, E, K) = 1. (2.12) 

Then as in [8] and [11], we have 

Pr{±oo, E) > 0, 



iPj{±oo, E) 

\ (-1)^+1 v/pij(±oo, E) 

fX poo 

/ PR{y,E)dy = {x - xi{E)) pr{oo, E) + / {pR{y, E) - pr{oo, E)) dy 

Jxi{E) Jxi{E) 

+ O {\x\~^~'^^ , as a; — )■ cxD, and 

PX P — OO 

/ PR(,y,E)dy = {x - xq{E)) pr{-oo, E) + / {pn^y, E) - pr{- 00, E)) dy 
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+ O [\x\ ^ as a; -oo. 

The error estimates here are uniform for E & A. 

We thus have an incoming wave asymptotically described for large negative x by 

!x(^E)iPR'yy^ E)-Pr(-co, E)) dy/h 

da (-00, E, h) ^=^= ^^^VR{-oo,E)x^{E)|h ^ipn{-oo,E)x/H_ 

^/Pr{-00, E) 

The reflected wave is asymptotically described for large negative x by 

p-* !x^E) (PRiV' E)-Pr{-co, E)) dy/h 

di{-oo, E, h) ^=^= g+iPfl(-oo,£)xo(£)/n g-*pB(-oo,i?)xA_ 

A/pi?(-oo, E) 

The transmitted wave is asymptotically described for large positive x by 

it^Epn'^y^ E)-Pr{+cx., E)) dy/h 



Q-~tPR[-oo,E)xi(E)/n ^ipR{+oo,E)x/n ^2 13) 

a/ph(oo, E) 

We obtain the solution ip{x^ E, h) that we use in (12.31) by normalizing the incoming 
flux. To do so, we simply divide the whole solution C{x, E, K) by the constant 

fx(TE)(PR(y' E)-pr{-oo, E)) dy/h 
d2(-00, E, h) -^=^= ^-ipR{-^.E)x,{E)/h_ 

/pi?(-oo, E) 



We obtain our main results by analyzing the tunneling wave by studying the large x and 
t asymptotics of ( 12. 3p with this '^{x^ E, h) and an any appropriate choice of Q{E, h). 

2.3 Complex WKB analysis 

We need to compute the asymptotic behavior, as a; — t- —00, of the solution to (12. 9p that 
satisfies (12.121) . We do this by applying the complex WKB method in order to avoid matching 
of the solutions at the real turning points xq{E) and xi{E), where the equation is ill-defined. 
So, we consider (12. 2p and (12. 9 p in the strip 5*0, in the complex plane containing the real 
axis, with possible branch cuts at xq{E) and xi{E). We now replace the variable x by 
z, to emphasize that the variable is no longer restricted to the real line. The solution to 
(12. 2p is analytic for z E Sa-, but the solution to (12.90 has singularities at the turning points. 
As Re 2; tends to +00 in S'q,, our assumptions on the behavior of the potential ensure that 
the asymptotic values of the coefficients Cj{z) are independent of Im^;. We can thus start 
the integration of (12.90 above the real axis at the extreme right of the strip, with asymptotic 
boundary data (12.120 . Also, our assumptions imply the existence of two paths from the 
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right end of the strip Sa to its left end, with one of them, 7^, passing above the two turning 
points, and the other, 7^, passing between them. We parameterize these with t E M. and 
assume they satisfy 

7#(t) € Sa, with hm Re 7#(t) = =Foo, and sup |7#('^^)| < 00, 

where # stands for a or b. Note that Im7a(t) > for all t, whereas Im7f, changes sign 
exactly once and is positive for t large and negative. Finally, and this is the main property of 
these paths, the imaginary part of J^_^(^^^ p{z', E) dz' along 7^ is decreasing. Such paths 
are called dissipative. 

The existence of these dissipative paths is proved as in [10] and [9]. Close enough to the 
real axis, there exist level lines of the function Im j^^f^^^ vi.^' 1 -^) dz' which are essentially 
parallel to the real axis for Rez > xi{E) and Re^ < xq{E). For xq{E) < Re^; < xi{E), 
these lines can be connected by means of level lines of Re v{z' 1 E) dz', along which 

Im /jf^j-^-) p{z', E) dz' is strictly decreasing. As a local analysis reveals, the connections can 
be made in a fashion without losing the dissipativity property. It is readily seen by 
inspection, that 7a can be constructed this way. For 7^, one starts as for 7^, and between 
xo{E) and xi{E), one uses a level line of Im f^^f^^^ p{z' , E) dz' to cross the real axis. Then 
one connects to a level line of Re p{z', E) dz' and proceeds as described above to 

connect to —00 below the real axis. 

We can integrate (12. 9p along these two different paths and compare the solutions for large 
negative Rez with z E Sa- These two integrations describe the same solution to (12. 2p since 
it is analytic. Moreover, the asymptotic values of the coefficients as /i — )■ can be controlled, 
because both these paths are dissipative. 

We choose two specific branches, p"'{z, E) and p^{z, E), of the multivalued function 
p{z, E). For p°'{z, E) we place vertical branch cuts below the real axis, extending down 
from Xq{E) and Xi{E). For p^{z, E) we place a vertical branch cut below the real axis 
extending down from Xi{E), and a vertical branch cut above the real axis extending up from 
Xo{E). These two functions are then uniquely determined in their respective regions Sa and 
S^ by the requirement that they both equal Pr{z, E) when z is real and greater than Xi{E). 
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The following are satisfied when z is on the real axis: 

^^(a;, -E) > 0, if X > xx{E) 

e'""/^ Pnix, E), if xo{E) < x < Xi{E) 

^ - pr{x, E) < 0, if X < xo{E). 



p%x, E) 



and 



/(x, E) 



Pr{x, E) > 0, if X > xi{E) 

e'""/^ Pr{x, E), if xo{E) <x < xi{E) 

Pr{x, E) > 0, if X < xo{E). 
For ^ = a, b, the function p'^lz, E) is analytic in a neighborhood of the path 7#. 

We define ^pfiz, E), to be the analytic continuation in 5* of the vector (pj{z, E) de- 
fined in (12. 7p . We also define / -E') c^y to be the analytic continuation in 5* of 

I p{y, E) dy, which is already specified for x > xi{E). 

Lemma 2.1 For any real z < xo{E), the following are satisfied: 

1 



I (p^{z) = (P2{Z) 



^Pr{z, E) 



xi{E) 



p\y, E)dy 



xi{E) 



Vpniz, E) 



xoiE) 

p\y,E)dy = i I pR{y,E)dy + I pR{y,E)dy 



xi{E) 



f 

Jxo{E) 



p\y, E)dy = -2i I pR{y, E)dy. 

xiiE) Jxo{E) 



xi{E) 



(2.14) 



Proof We simply follow the analytic continuations of p in the respective domains. 



Remarks 



rxi[a) 

i) The quantity 2 / E) dy can be expressed as a contour integral 

Jxn(E) 

p{z, E) dz > 0, 



pxiiE) 

K{E) = 2 / pR{y, E) dy 

JxoiE) 
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where 7 is a simple negatively oriented loop around the two turning points, and p{z, E) 
is the analytic continuation of Pr{x, E) for x > Xo{E). This shows that K{E) is 
analytic for E in a. complex neighbourhood of the energy window A. 

ii) Equation f l2.14p shows that the analytic continuations {c\{z), c\{z)) of the coefficients 
coincide with {di{z), d2{z)) in (12.111) for z < xq{E), up to multiphcative constants. 

Coming back to the differential equation (12. 9p . we denote the analytic continuations of its 
solutions in S** as {cf,cf). We consider the analytic function $ for z < xo{E) and the two 
different analytic continuations of its decomposition (12. 8p at z. These two representations 
must agree. This and Lemma [2.11 imply the following: 

Lemma 2.2 For any z < xo{E), we have 
i 4{z, E, h) e^(^)/^' = cl{z, E, h) 

I c\{z, E, h) e-^(^)/'' = c^(z, E, h). 

Remark The identities in the two previous lemmas are actually true for any z with 

Re 2; < Xq{E). 

The WKB analysis of (12.90 along the dissipative paths 7* and assumption (12. ip now 
yield the following lemma, as shown in (TU], [H], [H], [TT] - 

Lemma 2.3 We have 

cl{-oo,E,h) = 1 + OE{h), as h~^0, 

cl{—oo,E,h) = 1 + Oe{K), as /i — )■ 0, and 

cf (x, E, h) = cf(±oo, E, h) + Oe{1/\x\^^'') as x-^ ±00, 

where the remainder terms are analytic in E, for E in a complex neighborhood of the 

real set A. Moreover, —— C2(— 00, E, h) and the Oe{^/\x\^) are uniformly hounded 

dE 

as h ^ 0. 
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As a consequence of this lemma, for x ^ 1, we have 

P^(I^(E) iPRiV' E)-pRioo, E) dy)/h -ixi{E)pR{oo, E)/h 

C{x, E, h) = gip«(+oo,£;)x/s 

A/pij(oo, E) 

X (1+ Os(/i) + OE{\l{h\x\^^'')) 

and, for a; ^ — 1, we have 

pi<yrj^''\vR(.s.E)-VR{-oo,E)ds)lh ixo{E)pR{-oo,E)/h 
Cix, E,h) = -I e^(^)/(2r.) ^ ^ ^-^PRi-oo,E)./r^ 

VPij(-oo, E) 

p-i{f"t''\pRis>E)-pR{-oc,E)ds)/h -ixo{E)pR{-oo,E)/h 
^ ^K{E)/{2h) ^ ^. gipfl{-oo,i=;)x/?l 

^JpR{-O0, E) 

^K{E)/{2h) ^ipR{~oo,E)x/h ^ ^-ipRi-oc>,E)x/h 

X (Oe(/i) + OE{l/{h\x\'+'')). (2.15) 

2.4 Large Time Asymptotics of the Tunneling Wave Function 

We now consider the large time behavior of the transmitted wave packet. We denote this 
wave packet by t, h). We construct it as a time-dependent superposition of the 

normalized generalized wave functions ip{x, E, h), where x > 1 > max xi{E). 
The specific superposition we use is 

X{x, t,h)= [ Q{E, h) e-''^'^' ^{x, E, K) dE, (2.16) 

J A 

where for x > max x-i (E), 

EGA ^ ' 

ip{x, E, h) (2.17) 
g-x(E)/(2ro ^p(_oo, E) c^(x, E, h) e'^=^'i(^^) ^^^^'^^'^^/'^ 

^/p{x^ 4{-00, E, h) e-^/-2i'"(Pfl(2/.^)-P«(-°^.^))A (,-ipR(-oo,E)xo{E) ' 

See Remark ii) after Lemma 12.11 The asymptotics we have established show that for large 
positive x. 



p[+oo, E) 

X e-^'^(-^)/^' (1 + r(x, h)) dE, 
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where 

-xoiE) 



io{E) 



POO rXQiytj) 

/ (P niy, E) -p nioo, E))dy - / {p R{y, E) - pr{- oo, E)) dy 

Jxi{E) J -oo 



+ pr{-^,E){xi{E)-xo{E)). 
The error term r(x, E, K) in this expression satisfies 

r{x,E,h) = 0(h + + = 0(h+ ^ 



uniformly for E E A, x > 1 and h small enough. 

We prove below that xi^^ K) asymptotically propagates freely to the right for large 
positive t. 

For G A, we define 

a{E) = G{E) + K{E)/2 and k{E) = J{E) + uj{E), (2.18) 

where G and J are the functions in (12.41) . We then have 

X{x, t, h) (2.19) 



A V Py+^' ^! 

We obtain the small h asymptotics of this integral by Laplace's method. We first state a 
result concerning the large x and large t behavior of xi^-, t, h), whose proof follows from the 
methods of [8] and [11] , but is easier. The detailed analysis of the x and t dependence yields 
the following result. See [8] and [11] for details. 

Theorem 2.4 Assume the function a{E) has a unique non- degenerate minimum at 
E = E* m A. Define k{E) = pr{oo, E) and k* = k{E*). 

There exist 6 > and > 0, such that for t > Tfi and all x > 1, the transmitted wave 
satisfies 

X(x, t, h) = xcaussix, t, h) + 0{e-'^^^*^/'h'/'^') , 
where the error term is measured in the norm, uniformly for t > T;i, and 



XGaussix,t,h) = V2Vhk^ P{E\h) J ^l,^' Zl e""^""*^^' 



p(+oo, E*) 
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X 



exp{-i{tE* + k{E*) - k*x)/h} 



X exp 



(x - k*{t + k'{E*)))^ 
2h {^a{E{k))\k* +t{t + ^^{E{k))\k*)) 



Proof Outline The proof of this theorem is very technical, but is very similar to ones 
presented for Theorem 5.1 of [8] and Theorem 6 of |TT]. One computes the leading term 
Xcaussi^y t, h) by a rigorous version of Laplace's method, paying attention to the dependence 
of the remainder terms on the parameters x and t. The norm of this leading term is 

0' .3/. ^ ylg^)" ( f^l^aiEmY", (2.20) 

which is O ( h^^^ g-^C-^*)/'* ^ . By the methods of used in [8] and [IT], the norm of the error 
term induced by r{x, E, h) under the integral sign in f l2.19p is of order /i3/4+(5 
some (5 > 0, provided t is large enough. Note that Lemma 1 and Proposition 5 of [2J allow 
us to get better control of T^. (See below.) | 

2.5 The Transmitted Wave Function Shortly After Tunnehng 

Mimicking [2], we shall now address the behavior of the transmitted wave for finite values 
of X and t, shortly after the transmitted wave has left the region where it emerges from 
the potential barrier. Because the exponential decay of transmitted wave computed from 
the behavior of p{z,E) on the real axis appears as a factor, see (12.171) . the analysis of the 
semicassical behavior of the coefficient c"(x, E, K) for finite values of x > Xi{E) is simpler 
than in [2]. By contrast, the appearance for finite x's of the corresponding exponentially small 
factor for the above barrier refiection required to pass to the superadiabatic representation 
in P]. This is not necessary here so that we can stick to the adiabatic basis fl2.10p . We 
shall not, however, examine the more delicate details of how the transmitted wave actually 
emerges from the barrier. One should be able to address this much more technical topic by 
adapting the results of [2]. 
Let 

rxo{E) 

P{E) = - {pniy, E) - pr{- oo, E))dy - pr{-oo, E) xo{E), 

J —oo 
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S{x,t,E) = - [ PR{y,E)dy + p{E) + J{E) + Et, 

JxiiE) 



and 



In the region of moderately large positive x, but far from the potential barrier, the 
transmitted wave is described by the following theorem, which requires faster decay of the 
potential to its asymptotic value. 

Theorem 2.5 Let v > 21/2. There exist 5 > 0, t > 0, C > 0, and l3 > 0, such that 
for all t > T and sufficiently small h, 

{Xmodix, t, h) if I <X < Ch'^ 
X^aussi^,t,h) if Ch-^<x } + 0(;i3/^+^e--(^*)/'') 

otherwise 
where the error term is measured in the norm, 



2 

t.E*)) 

1 



^ ^ ^/a"{E*) + tS"{x, t, E*) 
and ' denotes the derivative with respect to E. 



Proof Outline We follow the main steps of the proof of the corresponding result for 
Theorem 5 of [2], with one notable exception. Since we do not use any superadiabatic 
representation, the next to leading order term in the asymptotics of c"(x, E, K) is of too low 
an order to be treated as in |2j. We briefly address this issue in more detail. By integration 
by parts, we see that for x > 1, 

c\x,E,h) = 1 + zh iS^l^dy + OEin'/x"^'). 

When we integrate against the energy density Q{E, h), the remainder term can be dealt 
with by using Lemma 1 of [2j. The non-zero next to leading term is only of order h/x^^^ 
and the error term it generates can be bounded as follows: 
Let ?7(x, t, h) = x{x, t, K) — t, h), where 

X{x, E, h) (2.21) 

E) 4(-00, h) e~'S-t''\pRiy^E)-PRi-^'E))/h ^-^p^i-oo,E)xo{E) ' 
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The error term whose norm we need to bound has the exphcit form 

r foo (dp/ 

g{x,t,h) := ih Q{E,h)e-^'^'/^x{x,E,h) / , 3/ 1 dy dE 

= h [ e-'^*/'^ Q{E, h) fix, E) e'^-'i(^) p^^^''^)'^^/^ dE, 

where /(x, E) = 0{l/x'^~^^), uniformly for E E A, Q{E, h) is independent of x, 
and \Q{E, h)\ behaves essentially as e'^^^*)!^ times a Gaussian in [E — E*)/\/h. (See 
fl2A8|) l We compute 

/ \g{x,t,h)\^ dx (2.22) 

Jx>l 

= QiE,h)Q{E',h) fix, E)J{x, E')e'^^ii^^^'"'^''^^-''''^''^'^^''''^^dxdEdE'. 

JaxA Jl 

Let < < 1. We split the integration range into the sets where \E — E'\ < and 
\E — E'\ > . We perform an integration by parts in x on the latter set to get 
fix, E) Jix, E') ^'L\,E)iPniy'^)-Pniy,E'))dy/n^^ 



i h 



- ipnil, E)+pnil, E')) /(I, E) /(I, E') e'A(^) ipniy,E)-p,(y^E'))ay/n 



2iE-E' 

+ 2iE-E') l ^{{PR{^.E)+pnix,E')) fix,E) fix,E')) 

X e*-^"i(-^) iPRiy,E)-PR{y,E'))dy/h 

The absolute values of both terms are bounded by Ch}~^ , where C is uniform in i?. So, they 
contribute a factor C/i^~^||Q||f to fl2.22p . Similarly, the integral in x in fl2.22p is bounded 
uniformly in E and, for some constant C, 

[ \QiE,h)\ \QiE',h)\ X{\E-E'\<no} dE dE' < C \\Q\\l 

JaxA 

Since A is compact, ||QI|i ^ l^l^'^/^IIQIh- Using this and the the estimate 
IIQII2 < Ce-"(■^*)/^/ll/^ we eventually get for 6 = 2/3, 

\\gh < C e-"(^*)/^^ ^1^/12 < C e-"(^*)/^' 

The rest of the proof, which consists of showing that x can be approximated by Xmod 
and Xcauss different values of x, now relies on Lemma 12.31 and on arguments identical to 
those in the proof of Theorem 5 of [2] . | 
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While explicit and concise, the approximation above does not make apparent where the 

transmitted wave is actually located. To have a better idea of the position of this wave 

d 

function, we define qt to be the unique solution in x to —— S{x, t, E*) = with qt > and 

oE 

qt > 0. The function qt is actually the classical trajectory in the potential V with energy 
E*, the velocity of which is bounded from above and below. We define 



f + k\ Po{qt, E*) t , -Of + T?*\\in\ 

Xa.U^,t,K) = ^„„(^.) ^ .,g„(^,_ , exp{-(a(g ) + ^S(., t, E ))/fi} 



{x - qtf 

X exp < - 



2hpR{qt, E*Y{a"{E*) + iS"{qt, t, E*)) 
This wave packet is a Gaussian that is centered on the trajectory qt, and whose width is of 
order ^/h. 

This leads immediately to the following corollary. (See Theorem 6 of [2].) 

Corollary 2.6 There exist Xq > and 6 > 0, such that for all times t, with 
Xq < qt < Ch^^ , we have, in the L"^ sense, 

X{x,t,h) = XGaussix, t, h) + 0(;i3/4+5 

where Wxcaussix, t, h)\\L2 = O ( ^^/^ e^^^*)/^ ) . 
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